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We conducted pre- and post-perihelion observations of Comet C/2009 PI (Garradd) on UT 201 1 October 13 
and 2012 January 8, at heliocentric distances of 1.83 and 1.57 AU, respectively, using the high-resolution 
infrared spectrometer (NIRSPEC) at the Keck II 10-m telescope on Mauna Kea, HI. Pre-perihelion, we 
obtained production rates for nine primary volatiles (native ices): H 2 0, CO, CH 3 OH, CH 4 , C 2 H 6 , HCN, C 2 H 2 , 
H 2 CO, and NH 3 . Post-perihelion, we obtained production rates for three of these (H 2 0, CH 4 , and HCN) 
and sensitive upper limits for three others (C 2 H 2 , H 2 CO, and NH 3 ). CO was enriched and C 2 H 2 was depleted, 
yet C 2 H 6 and CH 3 OH were close to their current mean values as measured in a dominant group of Oort cloud 
comets. This may indicate processing of its pre-cometary ices in a relatively oxygen-rich environment. 

Our measurements indicate consistent pre- and post-perihelion abundance ratios relative to H 2 0, sug- 
gesting we were measuring compositional homogeneity among measured species to the depths in the 
nucleus sampled. However, the overall gas production was lower post-perihelion despite its smaller helio- 
centric distance on January 8. This is qualitatively consistent with other studies of C/2009 PI , perhaps due to 
seasonal differences in the heating of one or more active regions on the nucleus. 

On October 13, the water profile showed a pronounced excess towards the Sun-facing hemisphere that 
was not seen in other molecules, including H 2 0 on January 8, nor in the dust continuum. Inter-comparison 
of profiles from October 13 permitted us to quantify contributions due to release of H 2 0 from the nucleus, 
and from its release in the coma. This resulted in the latter source contributing 25-30% of the total observed 
water within our slit, which covered roughly ±300 km by ±4500 km from the nucleus. We attribute this 
excess H 2 0, which peaked at a mean projected distance of 1300-1500 km from the nucleus, to release from 
water-rich, relatively pure icy grains. 

Published by Elsevier Inc. 


1. Introduction 

More than any other Solar System material, the ices contained 
in the nuclei of comets (i.e„ their native ices) * 1 retain a relatively 
well-preserved record of conditions when and where these primitive 
small bodies formed. Comets reside for long periods of time in the 
outer Solar System, and clues to its formation and evolution are pro- 
vided through their volatile compositions (e.g„ see reviews by Mum- 


* Corresponding author at: Goddard Center for Astrobiology, NASA Goddard 
Space Flight Center, Mail Stop 693.0, Greenbelt, MD 20771, USA. 

E-mail address: michael.a.disanti@nasa.gov (M.A. DiSanti). 

1 We use the term "native ices” to refer to those housed in the nucleus, which upon 
sublimation release primary volatiles into the coma. 

0019-1035/$ - see front matter Published by Elsevier Inc. 
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ma and Charnley, 2011; DiSanti and Mumma, 2008; Bockelee- 
Morvan et al., 2005). 

Comets today belong to (at least) two distinct dynamical reser- 
voirs, the Kuiper Belt (KB) and the Oort cloud (OC), thought to be 
the dominant source regions for ecliptic and nearly isotropic 
(long-period and dynamically-new) comets, respectively. Until rel- 
atively recently it was thought that short-period comets formed 
more or less in situ in the KB (at heliocentric distances R h > 30 AU), 
while those from the OC formed in the R h ~ 5-30 AU region and 
were subsequently ejected through gravitational interaction with 
the giant planets. If true, this suggested distinctly differing compo- 
sitions might be expected for these two dynamical populations of 
comets (Mumma et al., 1993). 

This reasoning has since evolved to new paradigms, stemming 
from the Nice model (Gomes et al., 2005; Morbidelli et al„ 2008) 
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Table 1 

Observing Log for C/2009 PI (Garradd). 3 


NIRSPEC setting 

UT start-end 

Molecule (order) b 

T int c (min) 

Pre-perihelion - UT 2011 October 13 
R h = 1.830 AU, A = 1.836 AU, dzl/dt = +19.2 km s _1 
KL2 

, dR/,/dt= -12.2 km s -1 
04:42-05:42 

H 2 0 (26) 

HCN, C 2 H 2 , NH 3 (25) 
CH 4 (23) 

H 2 CO (21) 

36 

MWA 

05:52-06:15 

H 2 0 (15, 16) 
CO (16) 

8 

KL1 

06:26-06:35 

H 2 0 (26, 27) 
C 2 H 6 (23) 
CH 3 OH (22, 23) 

6 

Post-perihelion - UT 2012 January 8 
R h = 1 .566 AU, A = 1 .852 AU, dzl/dt = -20.6 km s" 
KL2 

\ dRh/dt = 3.4 km s -1 
15:21:15:43 

H 2 0 (26) 

HCN, C 2 H 2 , NH 3 (25) 
CH 4 (23) 

H 2 CO (21) 

15 


3 On October 13 (January 8), flux calibration was established through spectra of the bright infrared standard star HR-8541 (HR-4357) acquired immediately following (prior 
to) the comet observations. 

b Measured primary molecules and NIRSPEC echelle order in which their emissions occur. 
c Total on-source integration times. 


that predicts considerable disruption of the outer Solar System due 
to outward migration of young Saturn through the 1:2 orbital res- 
onance with Jupiter. The more recent "Grand Tack" model (Walsh 
and Morbidelli, 2011) invokes inward migration of the giant plan- 
ets within the first 5 x 10 6 years of the Solar System - as a means 
of explaining the relatively small mass of Mars and also the forma- 
tion of small bodies (especially asteroids) in the inner Solar System 
(R h ~ 1-3 AU) - followed by their outward migration. These mod- 
els predict compositional “blurring" of comets along dynamical 
lines, in the sense that diversity in cometary compositions should 
exist in both KB and OC reservoirs, and indeed this is observed. 
The models however do not account for the observed discovery 
rate of nearly isotropic comets. As a result, it was argued that a 
large fraction (up to 90%) of OC comets may have originated in 
the environs of neighboring stars in the Sun’s birth cluster (Levison 
et al., 2010), although more recent modeling suggests revising this 
fraction downward (Brasser et al., 2012). 

Here we report abundances for nine primary volatiles in the 
dynamically new Comet C/2009 PI (Garradd; hereafter C/2009 
PI ) based on observations obtained near R h = 1.8 AU pre-perihelion 
and 1 .6 AU post-perihelion. We compare our results with other 
infrared studies of C/2009 PI, place them in the context of our 
ongoing compositional taxonomy among comets measured to date, 
and discuss possible implications regarding the processing history 
of the ices in this comet. We compare spatial emission profiles, and 
thereby estimate the fraction of H 2 0 released from what appear to 
be water-rich grains in the coma. 

2. Observations and data processing 

We observed C/2009 PI pre-perihelion on UT 2011 October 13, 
and also shortly following perihelion (December 23, q = 1.55 AU) 
on 2012 January 8 (Table 1 ) using the NIRSPEC infrared high-reso- 
lution (2 /Aa ~ 24, 000) spectrometer (McLean et al., 1998, 2000) at 
the Keck II 10-m telescope on Mauna Kea, Hawaii. On October 13, 
C/2009 PI reached the Keck II west elevation limit at 06:32 UT, 
which translated to approximately 1 h and 40 min of clock time 
following focus and acquisition at the start of the night. In spite 
of this limited available time, the combination of broad spectral 
coverage and high-throughput optics afforded by NIRSPEC allowed 
us to sample a large suite of primary volatiles using three distinct 
cross-disperser/echelle combinations (Table 1). On January 8, 


C/2009 PI became available in the east toward the end of the night, 
only allowing accommodation of a setting that encompasses the 
largest number of organic molecules simultaneously with H 2 0. 

For our observations, we oriented the slit along the Sun-comet 
line (PA = 76° on October 13, 330° on January 8), 2 which was main- 
tained throughout our observations of C/2009 PI by using NIRSPEC 
in active rotator (i.e„ ’Position Angle’) mode. Comet and supporting 
flux standard star spectral frames (Table 1, note ‘a’) were acquired 
in our standard ABBA cadence, with both A- and B-beams present 
in the slit (e.g., see Villanueva et al., 2011a; Bonev, 2005; DiSanti 
et al., 2006; Mumma et al., 2001a) and separated by 12 arcsec. 

Cometary and stellar data were processed using our custom algo- 
rithms (e.g., Villanueva et al., 2011a; Bonev, 2005; DiSanti et al., 
2001 ) that correct for anamorphic optics and provide highly accu- 
rate spatial and spectral registrations. We matched abundances of 
molecular absorbers in the terrestrial atmosphere to extracted spec- 
tra, using the Line-By-Line Radiative Transfer Model (LBLRTM; 
Clough et al., 2005) together with our custom updates (Villanueva 
et al., 2011b). The high terrestrial water burden on October 13 
(~5 precipitable mm) limited the number of accessible cometary 
water lines, particularly in order 27 of the KL1 setting (see Sec- 
tion 3.1). The water burden on January 8 was relatively modest 
(~1.5mm), however the air mass of C/2009 PI was higher com- 
pared with October 13. Spectral extracts, representing sums over 
nine spatial pixels (1.78 arcsec) centered on peak emission intensi- 
ties, are shown for October 13 (Figs. 1 and 2) and January 8 (Fig. 3). 

3. Rotational temperatures, production rates, and spatial 
distributions 

3.1. Excitation analysis 

Because vibrational bands are in general not sampled com- 
pletely, obtaining accurate production rates for observed mole- 
cules requires knowledge of the relative populations of their 
rotational levels. These are characterized by a rotational tempera- 
ture (T ro t), determined by comparing observed and modeled inten- 
sities (fluorescence efficiencies; i.e„ g-factors) for emission lines 


2 Position Angle (PA) is defined as the direction in the sky plane of the extended 
heliocentric radius vector measured eastward from north, termed “PsAng" in the JPL- 
Horizons system (http://ssd.jpl.nasa.gov/horizons.cgi). 
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Fig. 1. Spectral extracts for several echelle orders, with modeled continuum and 
molecular emissions in C/2009 PI on UT 2011 October 13. Residuals and ±1 it 
stochastic noise envelopes are also shown. 


spanning a sufficiently large range of rotational energies. Our line 
g-factors are based on ro-vibrational band fluorescence models 
we developed for each molecule: H 2 0 multiple bands (Figs, la, 
2c and 3c; Villanueva et al., 2012a), OH* Vt and v 2 - (Figs, lb, 
2c and 3a-3c; Bonev, 2005, Bonev et al., 2006), C 2 H 6 v 7 (Fig. 2a; 
Villanueva et al., 2011b), CH 4 v 3 (Figs, lb and 3b; Gibb et al., 
2003), C 2 H 2 v 3 and v 2 + v 4 + v 5 (Figs, lc and 3a), CO v 3 (Fig. la) 
(Villanueva et al., 2011a), NH 3 Vj and v 3 (Figs, lc and 3a; Villanueva 
et al., 2013; Kawakita and Mumma, 201 1 ), HCN v 3 (Figs, lc and 3a; 
Lippi et al., 2013; Villanueva et al„ 2013), H 2 CO Vj and v 5 (Reuter 
et al., 1989; DiSanti et al., 2006), and CH 3 OH v 2 (Fig. 2a; DiSanti 
et al., 2013) and v 3 (Fig. 2b; Villanueva et al., 2012b). These models 
make use of a realistic solar flux that includes Fraunhofer lines and 
incorporates the cometary heliocentric velocity (dR h /dt in Table 1) 
in calculating fluorescent pumping rates (see Villanueva et al., 
2011a for additional details). 

The high terrestrial water burden on October 13 resulted in 
poor transmittance for many of the cometary H 2 0 lines, and so 
introduced correspondingly large uncertainties in line-by-line pro- 
duction rates retrieved from transmittance-corrected line fluxes 


Fig. 2. Additional spectral extracts for C/2009 PI on UT 2011 October 13. The OH 
line indicated in panel c here and in Fig. 3c is a combination of fluorescent and 
prompt emission. 

(see Eq. (1) below). The H 2 0 line near 3514.4 cur 1 (rest frequency) 
is normally the strongest observed in the L-band, but it suffered 
from high opacity (transmittance less than 1%), and was not de- 
tected (Fig. 2c). Using observed lines in the three orders shown 
in Fig. 2c yielded a poorly constrained rotational temperature for 
water (T rot = 42 +26 /_ 24 I<). 

We retrieved more precise rotational temperatures independently 
for HCN, CH 4 , CO, C 2 H 6 , and CH 3 OH on October 13 (Table 2). Their 
weighted mean value (T rot = 48±4K; lcr uncertainty) is consistent 
with r rot measured for C/2009 PI near R h = 2.1 AU (40 ± 7 K; Villanueva 
et al., 2012c) and 2.0 AU (50 K; Paganini et al., 2012). We take our 
weighted mean Trot (48 K) to be the most probable temperature, 
and adopt this for all molecules on October 13. From our January 
8 observations, we measured Trot = 55 ± 10 I< based on water lines 
in order 26 (Fig. 3c). We were unable to determine T rot for HCN 
nor CH 4 , and therefore adopt 55 K for all molecules measured on 
January 8. 

We note that individual values of Trot measured on October 1 3 
agree within their lcr uncertainties (Table 2), as is generally found 
among molecules from studies of previous comets. This is 
consistent with electrons and fast H-atoms, released primarily 
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coma contained in the beam assuming release solely from the nucleus 
(i.e., purely native release). It is a geometrical parameter that de- 
creases with increasing distance from the nucleus and depends on 
R h , A, and gas outflow speed v g3tS , assumed to be spherically-symmetric 
and equal to 800 R,/ 05 m s -1 (Biver et al„ 2006, 3 2002). Because /(x) 
varies approximately as ty 1 for our field of view, our production rates 
are relatively insensitive to assumed photo-dissociation lifetime 
(e.g., see DiSanti et al., 2001). In practice Qn C for a molecule is the 
mean of its line-by-line nucleus-centered production rates weighted 
by their individual stochastic errors. 

For each molecule having emission with sufficient signal-to- 
noise, we then calculate the production rate at regular intervals 
(steps) along the slit (i.e., along the spatial profile) using Eq. (1), 
based on the encompassed Fi ine and incorporating the value of 
f[x) at each step. Outside the seeing-dominated inner coma, the 
production rates so determined increase to a constant value (with- 
in noise) that we take to be the terminal (or total) production rate 
(Qtot). This establishes a “growth factor” (GF) as the ratio Qtot/Qnc 
(e.g., see Appendix B2 in Bonev et al., 2006, and Dello Russo 
et al., 1998 for details regarding the “Q-curve” formalism), or: 

Q-tot = Qnc * GF. ( 2 ) 



Fig. 3. Same as Figs. 1 and 2, but for January 8. 


through photo-dissociation of H 2 0, controlling T rot in the coma (see 
Xie and Mumma, 1996; DiSanti et al., 2001 and references therein). 
It may also explain the similar T rot (also within ler) measured post- 
perihelion for H 2 0, for example if the smaller R h on January 8 was 
compensated by the lower water production rate Q(H 2 0). Our low- 
er post-perihelion Q(H 2 0) is consistent with pre-/post-perihelion 
asymmetries found by other observers of C/2009 PI (Section 4.3), 
and holds regardless of the exact values we adopt for T rot , both 
on October 13 and January 8 (see Section 3.2.1). 

3.2. Production rates 

Our formalism for determining production rates is based on the 
spatial distribution of emissions along the slit (i.e., the spatial pro- 
files; Fig. 4) and is well documented in the literature (e.g., Villanu- 
eva et al., 2011a; Bonev, 2005). A brief summary is provided here. 

We calculate a nucleus-centered molecular production rate 
(One) from signal contained in the central nine spatial pixels cen- 
tered on the peak emission intensity. For a given emission line 
(or, in the case of blends, for each group of lines); 


We generally find growth factors among simultaneously-mea- 
sured primary volatiles to agree within uncertainties. For this reason 
we adopt a common value, given by the weighted mean of measured 
growth factors for all (native) species within each given setting (see 
Section 3.2.1 ). However, release of molecules in the coma can result 
in larger growth factors. In Section 3.3 we discuss the influence of an 
additional, extended source of H 2 0 on October 13. 

3.2.1. Dependence of production rates on T rot and GF 

We obtained consistent production rates whether we used a 
common rotational temperature or each individually measured Trot 
(see note ‘b’ of Table 2). Including aT rot increases uncertainties in 
Qnc (listed in parentheses), and using individually measured aT rot 
likely overestimated crQnc since their T rot were less well con- 
strained. For this reason, and because Trot was not measured for 
all species, we take erQn C (and aQt ot ) based on the “combined” oT rot 
(±4I< on October 13) to provide the most consistent measure of 
molecular production rates. 

The influence of growth factor is also summarized in Table 2. 
We found no evidence that any molecules other than H 2 0 on Octo- 
ber 13 were subject to effects other than atmospheric seeing. Since 
seeing-related flux loss is a common systematic effect within a gi- 
ven setting, we used the same GF for all co-measured species ex- 
cept H 2 0 on October 13. For KL1 we applied the weighted mean 
GF from C 2 H 6 and CH 3 OH, emphasizing that their growth factors 
agree well within error. Similarly, for KL2 we used the weighted 
mean GF from HCN and CH 4 on October 13, and that from HCN, 
CH 4 , and H 2 0 on January 8. On the other hand, in addition to seeing 
the spatial distribution for H 2 0 on October 13 was seemingly influ- 
enced by release from a source in the coma, as evidenced by its lar- 
ger growth factors. In summary, for all molecules, both their 
production rates and abundances relative to H 2 0 were insensitive 
to individually measured rotational temperatures or (excepting 
H 2 0 on October 13) growth factors. 


_ 4nA 2 Fiine 
Tl /(*) g 1 


In Eq. (1), t 4 (s) and g 4 (W molecule -1 ) represent the photo-dissocia- 
tion lifetime and fluorescence g-factor (both evaluated at R h = 1 AU), 
F^ne (W m -2 ) is the line flux at the top of the terrestrial atmosphere, 
and geocentric distance (A) is expressed in meters. For each molecular 
species, the quantity f[x) represents the fraction of all molecules in the 


3.3. Significance of spatial profiles in assessing volatile release 

Comparing spatial profiles permits testing the nature of release 
of individual volatile species and dust. Fig. 4 shows emission pro- 
files from October 13 and January 8 for each molecular species 


3 The best-fit heliocentric dependence of % as for the 4 comets in Biver et al. (2006) 
is R h -°- 53±0 - 12 km s” 1 . 
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Table 2 

Volatile composition of C/2009 PI. 


Molec. 

Trod (K) 

Qnc b (10 25 s-’) 

GF C 

Qtot d (10 25 s ’) 

X/H 2 0 (%) e 

Pre-perihelion - 

UT 201 1 October 13 





KL2 setting 






h 2 o 

(48 ± 4) 

4640 ±312(352) 

1.94 ±0.069 

8980 ±752 

100 

HCN 

48 ±4 

14.4 ±0.84(0.95) 

1.498 ±0.071 

21.5 ±1.76 

0.23610.018 


55 +6 /_ 7 

13.3 ±0.99(1.45) 


20.0 ± 2.37 

0.219 ±0.0307 


55 +6 /_ 7 

13.3 ±0.99(1.45) 

1.497 ±0.086 

20.0 ± 2.42 

0.219 ±0.0312 

ch 4 

48 ±4 

52.9 ± 2.40(6.07) 

1.498 ±0.071 

79.2 ± 9.84 

0.882 ± 0.089 


50 ±10 

55.7 ±2.53(14.3) 


83.4 ±21.8 

0.913 ±0.249 


50 ±10 

55.7 ±2.53(14.3) 

1.501 ±0.148 

83.6 ±23.0 

0.914 ±0.261 

h 2 co 

(48 ± 4) 

6.53 ± 2.37(2.38) 

(1.498 ±0.071) 

9.78 ±3.60 

0.11 ±0.041 

nh 3 

(48 ± 4) 

29.7 ± 9.30(9.56) 

(1.498 ±0.071) 

44.5 ±14.5 

0.5010.16 

c 2 h 2 

(48 ± 4) 

3.36 ±1.32(1.32) 

(1.498 ±0.071) 

5.03 ±1.99 

0.056 1 0.023 

MwideA setting 






H 2 0 

(48 ± 4) 

5200 ±1680(1690) 

1.83 ±0.126 

9540 ±3150 

100 

CO 

48 ±4 

565 ± 26.6(35.9) 

1.52 ±0.056 

860 ± 63.0 

9.1210.80 


48 +10 /_9 

565 ± 26.6(63.2) 


860 ±101 

9.12 ± 1.22 f 

KL1 setting 






H 2 0 

(48 ± 4) 

6250 ±566(579) 

1.680 ±0.085 

10,600 ± 1100 

100 

c 2 h 6 

48 ±4 

45.9 ± 2.43(3.79) 

1.482 ±0.085 

68.1 ±5.61 

0.644 1 0.098 


44 +9 /_ 7 

44.7 ± 2.32(4.29) 


66.2 ± 6.36 

0.621 ±0.100 


44 +9 /_ 7 

44.7 ± 2.32(4.29) 

1.480 ±0.090 

66.1 ± 6.48 

0.626 ±0.102 

ch 3 oh 

48 ±4 

153 ±15.6(15.7) 

1.482 ±0.085 

226 ± 26.6 a 

2.14 1 0.38 


45 +6 /_ 5 

151 ± 15.5(15.7) 


224 ± 26.3 

2.12 ±0.38 


45 +6 /_ 5 

151 ± 15.5(15.7) 

1.501 ±0.266 

232 ±51.1 

2.15 ±0.52 

Post-perihelion - 

UT 2012 January 8 





KL2 setting 






H 2 0 

55 ±10 

4907 ±185(546) 

1.478 ±0.131 

7253 ± 698 

100 


55 ±10 

4907 ±185(546) 

1.459 ±0.178 

7159 ± 1183 

100 

HCN 

(55 ±10) 

9.2 ±0.90(1.83) 

1.478 ±0.131 

13.6 ±1.79 

0.1910.020 


(55 ±10) 

9.2 ±0.90(1.83) 

1.569 ±0.269 

14.4 ±3.79 

0.20 ± 0.063 

ch 4 

(55 ±10) 

35.7 ±7.35(11.5) 

1.478 ±0.131 

52.8 ±11.8 

0.7310.152 


(55 ±10) 

35.7 ±7.35(11.5) 

1.428 ±0.277 

51.0 ±19.2 

0.71 ± 0.293 

h 2 co 

(55 ±10) 

<4.9 (3o-) 

(1.478 ±0.131) 

<7.2 

<0.10 (3(7) 

nh 3 

(55 ±10) 

<35.7 (3(7) 

(1.478 ±0.131) 

<52.8 

<0.73 (3(7) 

c 2 h 2 

(55 ±10) 

<3.8 (3(7) 

(1.478 ±0.131) 

<5.6 

<0.077 (3(7) 


a Relatively precise rotational temperatures were measured for five primary volatiles on October 13 (HCN, CH 4 , CO, C 2 H 6 , and CH 3 OH). Their weighted mean (48 ±4 K, 
adopted for all molecules), and (in a second row) individually measured Trot are listed. On January 08, the value measured for H 2 0 was adopted for all molecules. Here and in 
subsequent tables, all uncertainties represent lcr and upper limits represent 3a. Values in parentheses indicate that temperatures were adopted from measurements of other 
molecules. 

b Nucleus centered production rate. For each entry, the first uncertainty corresponds to the larger of the stochastic or standard error among lines sampled, and the second 
entry (in parentheses) incorporates uncertainty in Trot- Production rates for individually measured Trot are also listed (second row). 

c Scaling factor (growth factor, GF) for Q outside the seeing-dominated region, relative to Qn C . Values in parentheses indicate that GF was adopted from measurements of 
other molecules. 

d Total production rate (Eq. (2)). For each molecule, the uncertainty includes both those in mean GF and T rot , with GF from October 13 based on HCN and CH 4 in KL2, and on 
C 2 H 6 and CH 3 OH in KL1. For these four molecules, the third row of entries includes uncertainties in both individually measured Trot and GF - in all cases Qn C and CW on 
October 13 agree (within 1 cr uncertainty) with their values at 48 I< and (for Qtot) with their values using the relevant mean growth factors. 

e Abundance, corresponding to the ratio of molecular production rates Qtot/Qtot(H 2 0). All molecules are assumed to have the same Trot. hence both pre- and post-perihelion, 
aT rot is not included in determining our “most probable” abundances (shown in bold). On October 13, uncertainties in GF are included since emissions span multiple settings. 

f CO abundance, compared to the weighted mean global production rate for H 2 0 among the three settings on October 13 (=9430 ± 606 x 10 25 molecules s '), at 
Trot = 48 ± 4 K (first row) and at the T ro , measured for CO. 

g The production rate listed for CH 3 OH is the weighted mean of those from orders 22 (232 ± 32 x 10 25 molecules s 1 ) and 23 (216 ± 39 x 10 25 molecules s '). 


(heavy trace in each panel), the dust continuum profile (light trace, 
co-measured within each order), and a representative point- 
spread-function (PSF, dotted trace) based on the measured stellar 
profile. In each order, the continuum profile defines the zero posi- 
tion for projected distance along the slit, and although distinct 
from co-measured molecular emission profiles as is generally the 
case for comets, it ensures uniform spatial registration for all tar- 
geted species. It also provides a reference position for extracting 
spectra and establishing nucleus-centered line intensities for spe- 
cies that lack a meaningful spatial profile owing to insufficient sig- 
nal-to-noise, i.e„ for H 2 CO, C 2 H 2 , and NH 3 , both on October 13 and 
January 8. 

The spatial profiles for H 2 0 on October 13 were clearly broader 
than those of other molecules, hence their larger growth factors. 
This suggests the presence of two ice phases, polar (i.e„ water-rich) 
and a-polar (CO-rich), as was discussed previously based on high- 
resolution IR spectra of C/2009 PI (Paganini et al„ 2012; Villanueva 
et al., 201 2c; see Section 4.4, and Appendix A). We note that the GF 


for H 2 0 in KL2 was larger than that in KL1 (compare water profiles 
in Figs. 4a and 4e). This could indicate a change in the projected 
(line-of-sight) position of H 2 0-rich (polar) grains in the coma and 
included within our (narrow) slit. It may or may not be associated 
with the reported nucleus rotation period of 10.4 h based on post- 
perihelion observations of C/2009 PI with D1XI/MRI on three UT 
dates spanning R h = 1.7-2. 1 AU (Bodewits et al., 2012). 4 

For H 2 0 in the MWA setting, we adopt the mean growth factor 
from KL1 and KL2, since we could not obtain a meaningful spatial 
profile due to insufficient signal-to-noise for water lines in this set- 
ting (see Fig. la). The MWA observations were conducted (tempo- 
rally) between KL2 and KL1 (Table 1), so this approach should be 
reasonably representative of any evolution in water release (as 
projected in our slit) over the course of our October 13 observa- 


4 If in a “complex” rotation state, the actual period for C/2009 PI would be a 
multiple of the reported 10-h period. 
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Fig. 4. Pre- (panels a-g) and post-perihelion (panels h-j) emission profiles in C/2009 PI, showing phase angle (/f) and slit orientation in the schematics at top. Heavy solid 
traces: Molecular emissions with ±1 a stochastic noise (dashed lines). Light solid traces: Continuum - in each panel, the centroid position of its co-measured profile defines the 
zero-point distance. Light dotted traces: PSF. Growth factors (GF) are listed for each molecular profile. The pronounced asymmetry displayed by Fi 2 0 on October 13 is not seen 
on January 8. 


tions. The value of Qt ot for H 2 0 in MWA is consistent with those in 
KL1 and KL2, albeit with larger uncertainty. 

C/2009 PI is the third recent comet to show evidence for (at 
least) two distinct phases of ice within the nucleus, one dominated 
by polar ice and the other by a-polar ice. In C/2007 W1 (Boattini), 
the polar molecules H 2 0 and CH 3 OH showed asymmetric profiles 
in the anti-sunward hemisphere, and a source of small particles 
composed of relatively pure polar (water-rich) ice was proposed 
(Villanueva et al., 2011a). Comet 103P/Hartley-2 showed symmet- 
ric H 2 0 and CH 3 OH profiles, while C 2 H 6 was strongly asymmetric 
and HCN was intermediate (Mumma et al., 201 1 ). Although a polar 
molecule the distribution of HCN was closer to that of C 2 H 6 in both 
C/2007 W1 and 103P, suggesting a source associated primarily 
with the a-polar ice contained in the nucleus. This also appears 
to be the case with C/2009 PI. 

In contrast, on January 8 the water profile showed no such evi- 
dence of a bimodal distribution, and its growth factor was consis- 
tent with those measured for CH 4 and HCN (Fig. 4h-j) and with 
release solely from the nucleus. For this reason we adopted the 
weighted mean GF for H 2 0, CH 4 , and HCN (1.48 ±0.13; Table 2) 
for all co-measured volatiles on January 8. Both H 2 0 and CH 4 
peaked slightly anti-sunward of the continuum, however their 
noise levels make it difficult to determine if any displacement is 
significant. We present additional pre-/post-perihelion compari- 
sons in Section 4.3. 


4. Comparisons among C/2009 PI observations 

4.1. Compositional comparison with other infrared observations 

Fig. 5 shows a comparison of our abundance ratios with those ob- 
tained near 2.1 and 2.0 AU (Villanueva et al., 2012c; Paganini et al., 
2012, respectively), relative to H 2 0 (Fig. 5a; Table 2) and HCN 
(Fig. 5b; Table 3). Those studies found growth factors of 1.8 for 
H 2 0 and 1.4 for other species, similar to those we found here. Our 
abundances on October 13 (at R h = 1.83 AU pre-perihelion) are con- 
sistent with those reported approximately 1 month earlier (near 
R h = 2A AU; Villanueva et al., 2012c) for HCN, CH 4 , NH 3 , and C 2 H 6 
(using NIRSPEC) and for CO (using CSHELL), whether expressed rela- 
tive to H 2 0 or HCN. Our retrieved abundances relative to H 2 0 are sys- 
tematically lower than those reported for CO, CH 3 OH, CH 4 , C 2 H 6 , and 
HCN from CR1RES observations near 2.0 AU, however it was pointed 
out (Paganini et al., 2012) that adopting the value for Q(H 2 0) from 
optical observations of OH (Section 4.2) brings the CRIRES results into 
the "normal” range (see Section 5). Perhaps the outgassing pattern 
(particularly of H 2 0) was changing with time, as suggested by 
comparisons among spatial profiles from infrared observations (e.g., 
compare Fig. 4a and e, and see Section 4.4). This could lead to abun- 
dance ratios relative to H 2 0 that are sensitive to rotational phase and 
slit geometry. Our measured abundances relative to HCN (Fig. 5b; 
Table 3) show much better agreement with the CRIRES results. 
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Heliocentric Distance (AU) DiSanti_etal_Fig5 

Fig. 5. Abundances of several primary volatiles in C/2009 PI (Garradd), expressed 
in percent relative to H 2 0 (panel a) and HCN (panel b) and retrieved using high- 
resolution IR spectroscopy. Note that for clarity, the points for CO and C 2 H 6 are 
offset by a factor of 10 relative to their measured values. Error bars denote ±la and 
upper limits (denoted by downward-facing arrows) are 3a. Perihelion 
(q = 1.5505 AU) is indicated by vertical dashed lines, and negative heliocentric 
distances denote pre-perihelion dates. Points (or groups of points) are labeled with 
instrument used: NIR = NIRSPEC/Keck II, CSH = CSHELL/1RTF, CRI - CRIRES/VLT. 
Measurements near R h = -2.1 AU are from Villanueva et al. (2012c), and those 
near -2.4 and -2.0 AU are from Paganini et al. (2012). The lower limit for the 
abundance of C 2 H 6 at -2.4 AU corresponds to the upper limit (3<r) found for H 2 0, 
which had not yet fully activated at that heliocentric distance. 

4.2. Comparison with optical observations of C/2009 PI 

D. Schleicher (private communication, 2011) reported water pro- 
duction rates (1.5 x 10 29 molecules s _1 ) from narrow-band photom- 
etry at R h = 1.97 AU (September 24) and 1.80 AU (October 18 and 20) 
that are approximately 50% larger than the values reported from the 
IR studies (Paganini et al., 2012; Villanueva et al., 2012c, and our Qtot 
for H 2 0 on October 13). This suggests release of additional water not 
encompassed by the (very narrow) IR slits. 

4.3. Pre-/post-perihelion comparisons 

Six of nine primary volatiles measured in C/2009 PI on October 
13 were also measured on January 8. Abundance ratios of HCN and 


CH 4 (relative to H 2 0) were unchanged (within 1 a uncertainties), 
and upper limits (3er) for H 2 CO, NH 3 , and C 2 H 2 on January 8 were 
also consistent with their measured values on October 13. Come- 
tary nuclei are thought to be processed to depths of tens of centi- 
meters during their residence time in the Oort cloud, primarily 
through irradiation by interstellar cosmic rays (Stern, 2003; also 
see Mumma et al., 1993). A typical comet loses approximately 
1 m from its surface per perihelion passage - for dynamically 
new comets this could lead to pre-/post-perihelion compositional 
differences resulting from surface erosion of irradiated layers. 
Regardless of the dynamical status of C/2009 PI (dynamically 
new or long-period returning), the agreement between pre- and 
post-perihelion abundance ratios suggests, for the active region) s) 
and to the depths sampled, compositional homogeneity among 
measured ices that may reflect its primordial endowment. 

A notable result is that the overall gas production in C/2009 PI 
was lower on January 8 despite being closer to perihelion: Qtot for 
H 2 0 was lower by ~23% relative to its mean value among settings 
on October 13 (Table 2, see note ‘f). Subsequent observations of 
C/2009 PI by EPOXI/MRI on UT 2012 February 23 (R,, = 1.7AU) 
and April 06 (R h = 2.1 AU) yielded values for Q(OH) (2.3 and 
1.4 x 10 28 molecules s _1 , respectively; Bodewits et al., 2012) that 
indicated a pronounced post-perihelion decrease in activity, at 
least in water production. SOHO/SWAN observations showed a 
steep (R/ 32 ) post-perihelion decrease in Q(H 2 0). This followed 
larger, relatively constant pre-perihelion production of H 2 0 from 
Ri, ss 2.1 to 1.6 AU (Combi et al, 2012) that was approximately 
three times our measured water production on October 13, 
consistent with additional release of H 2 0 in the much larger 
SOHO field-of-view and again suggesting aperture-dependent 
pre-perihelion water production. The lack of a “water excess” in 
our January 8 observations (compare Fig. 4h with Figs. 4a and 4e, 
and see Section 4.4) suggests loss of one or more active regions, 
perhaps due to seasonal changes in pole orientation affecting 
insolation received locally on the nucleus. The agreement among 
our pre- and post-perihelion abundance ratios (Table 2) is also 
consistent with loss of one or more sources releasing volatiles 
directly from the nucleus along with extended release of H 2 0, such 
that the production rates of H 2 0 and measured trace species 
decreased by similar amounts. 

4.4. Comparison of spatial profiles - establishing pre-perihelion 
"excess" H 2 0 emission 

On October 13, the weighted mean intensities for both H 2 0 and 
CH 3 OH were in the sunward-facing hemisphere (Fig. 4), yet only 
H 2 0 revealed a pronounced excess intensity that peaked between 
1300 and 1500 km sunward as projected on the sky plane. A weak- 
er secondary peak near +3000 to +3500 km was also seen in the 
KL2 water profile. This approximate “distance doubling" of peaks 
in the water profile could be associated with the same grain source, 
for example if material released previously and tied to nucleus 
rotation was also responsible for the weaker peak. Alternatively 
this secondary peak could indicate an additional source of 


Table 3 

Abundance ratios relative to HCN in C/2009 PI (Garradd) from IR spectroscopy. 3 


Reference (R h ) b 

CO 

c 2 H 6 

ch 4 

CH 3 OH 

c 2 h 2 

Villanueva et al. (2.1) 

37.1 ±6.7 

<3.14 

2.93 ± 0.68 


<0.36 

Paganini et al. (2.0) 

34.2 ± 3.8 

2.81 ±0.32 

3.56 ±1.78 

10.8 ± 1.73 

<0.36 

This work (1.8) 

38.0 ±4.6 

2.67 ± 0.47 

3.67 ± 0.48 

9.92 ±1.57 

0.23 ±0.10 

This work (1.6) 



3.64 ±1.03 


<0.39 


3 Entries represent 100 x Q(X)/Q(HCN) as shown in Fig. 5b, with uncertainties denoting la and upper limits denoting 3a. 

b Reference publication, and approximate heliocentric distance (AU, in parentheses). The first 3 rows in the table are pre-perihelion measurements, and the fourth row is 
post-perihelion. 
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Table 4 

Abundance ratios for reactants and products of H-atorn addition on interstellar grains. 


No." 

Comer 3 

Origin b 

CO/H 2 O a 

C 2 H 2 /H 2 O a 

CH 3 OH/H 2 O a + 

C 2 H 6 /H 2 O a 

1 

C/2001 A2 (LINEAR) 

NIC/OC 

3.9 ± l.l d 

0.5 ± 0.1 d 

3.0 ± 0.1 d 

1.7 ± 0.2 d 

2 

C/1996 B2 (Hyakutake) 

NIC/OC 

14.9 ± 1.9 e 

0.16 ±0.08 

1.7 ± 0.4 g 

0.62 ± 0.07 h 

3 

C/1995 01 (Hale-Bopp) 

NIC/OC 

12.0 ±0.4' 

0.31 ±0.10> 

1.7 ± 0.4 k 

0.56 ± 0.04' 

4 

C/1999 HI (Lee) 

NIC/OC 

1.8 ± 0.2 m 

0.27 ± 0.03 m 

2.02 ± 0.48 111 

0.67 ± 0.07 m 

5 

1 53P/Ikeya-Zhang 

HTC/OC 

4.7 ± 0.8 n 

0.18 ±0.05° 

1.75 ± 0.35 n 

0.62 ± 0.1 8 P 

6 

1 P/Halley 

HTC/OC 

3.5 q 

0.3 q 

1.7 ±0.4 

0.4 q 

7 

C/2004 Q2(Machholz) 

NIC/OC 

5.1 ± 0.51 s 

0.09 ± 0.01 s 

1.58 ± 0.06 s 

0.56 ± 0.03 s 

8 

8P/Tuttle H 

HTC/OC 

0.45 ± 0.09 

<0.04 u 

1.94 ± 0.1 0 U 

0.25 ± 0.03 t u 

9 

C/2001 WM, (LINEAR) 

NIC/OC 

0.52 ± 0.1 2 V 

<0.05 v 

1.12 ± 0.07 v 

0.47 ± 0.03 v 

10 

C/1999 S4 (LINEAR) 

NIC/OC 

0.9 ± 0.3 W 

<0.1 2 W 

<0.1 8 W 

0.11 ±0.02 w 

11 

73P/S-W 3-C H 

JFC/KB 

0.5 ± 0.13 x 

0.23 ± 0.06 v , <0.033 z 

0.14 ± 0.02 z 

0.18 ± 0.03 y,z 

12 

1 03P/Hartley 2 H 

JFC/KB 

0.15-0.45 A 

0.110±0.011 BC 

1.78 ± 0.08 b,c 

0.73 ± 0.03 B-C 

13 

C/2007 N1 (Lulin) 

NIC/OC 

2.19 ±0.10° 

0.076 ± 0.008° 

3.6 ±0.18° 

0.65 ± 0.08° 

14 

C/1999 Tl(McNt-Hartley) 

NIC/OC 

17 ± 4 C 

0.13 ±0.03- 

3.0 ±0.3° 

0.67 ± 0.033 e 

15 

C/2007 W3 (Boattini) 

NIC/OC 

4.5 ± 0.51 f 

0.29 ± 0.02 f 

3.67 ±0.11 

1.96 ± 0.04 1 

16 

C/2009 PI (Garradd) 

NIC/OC 

9.12 ±0.80'' 

0.056 ± 0.023 g 

2.1 ± 0.4 g 

0.64 ± 0.098° 


a ID numbers for comets refer to points in Fig. 6. Abundances are relative to F1 2 0 = 100. 

b Present dynamical classification (NIC = Nearly Isotropic Comet, HTC = Halley Type Comet, JFC = Jupiter Family Ecliptic Comet) and likely storage reservoir (OC = Oort 
cloud, KB = Kuiper Belt). 

c Previously published values for QfCH jOH) from IR observations using the v 3 Q-branch intensity have been inter-normalized to g(Q-br) = 1.5 x 10" 5 ph s" 1 molecule' 1 
(Villanueva et al„ 2012b). 
d Magee-Sauer et al. (2008b). 
e DiSanti et al. (2003). 
f Magee-Sauer et al. (2002a). 
g Biver et al. (1999a). 
h Dello Russo et al. (2002a). 

1 DiSanti et al. (2001). 

1 Magee-Sauer et al. (1999). 
k Biver et al. (1999b). 
m Mumma et al. (2001b). 
n DiSanti et al. (2002). 

° Magee-Sauer et al. (2002b). 
p Dello Russo et al. (2002b). 
q Eberhardt (1999). 

’ Eberhardt et al. (1994). 

! Bonev et al. (2009). 

' Bohnhardt et al., (2008) 

“ Bonev et al. (2008). 
v Radeva et al. (2010). 
w Mumma et al. (2001a). 
x DiSanti et al. (2007). 
y Villanueva et al. (2006). 
z Dello Russo et al. (2007). 

A Weaver et al. (201 1 ). 

B Dello Russo et al. (2011). 
c Mumma et al. (2011). 

D Gibb et al. (2012). 

E Mumma et al. (2001c) and DiSanti et al. (in preparation); 

F Villanueva et al. (2011a). 

G This work. 

H For 8P, 73P, and 103P, entries having multiple references represent weighted mean values among cited papers. 


water-rich grains having a larger characteristic scale for release of 
H 2 0 in the sunward-facing hemisphere. 

In contrast, CH 4 , CO, C 2 H 6 , and HCN exhibited similar distribu- 
tions, albeit with subtle differences: CH 4 and CO peaked (slightly) 
toward the anti-sunward-facing hemisphere and showed slight 
degrees of asymmetry in that direction, perhaps indicating 
“night-side” release. C 2 H 6 and HCN peaked coincident with the 
continuum and overall were more symmetric, however a statisti- 
cally-significant excess in HCN emission was seen near +3000 km 
(Fig. 4b), close in position to the simultaneously-measured second- 
ary water peak in the KL2 setting (Fig. 4a). Our data do not permit 
discerning if these two features are related. Addressing this and 
related questions requires detailed knowledge of the nature 
and actual three dimensional distributions of sources, outflow, 
and release. 

We estimated the fractional abundance of H 2 0 attributable to 
the “excess” emission through comparison with the combined 


profiles of C 2 H 6 and HCN. We postulate its origin in H 2 0-rich 
particles in the coma that contributed 25-30% of the total water 
within our slit (see Appendix A for details). 

Evidence for excess H 2 0 was also observed in C/2009 PI using 
N1RSPEC on UT 2011 September 9 (p = 26°, R h = 2.1 AU; Villanueva 
et al., 2012c). For those observations, the slit was oriented approx- 
imately 45° from the Sun-comet direction. The excess was most 
pronounced in the sunward-facing hemisphere (see Fig. 2C in 
Villanueva et al., 2012c) with a peak near projected distance 
+1500 km. A relatively smaller excess was seen on the opposite 
hemisphere, qualitatively consistent with our observations given 
the difference in slit orientations. A pronounced sun-side asymme- 
try was also seen in spectra of C/2009 PI near R h = 2.0 AU, using 
CRIRES with the slit along the projected Sun-comet direction 
(Paganini et al., 2012). Taken together, this suggests that the source 
region(s) on the nucleus giving rise to this excess H 2 0 may have 
persisted and was/were illuminated similarly for at least the 
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Fig. 6. Comparison of abundance ratios among comets (denoted by number in 
Table 4), related to processing of reactant (CO, C 2 H 2 ), and product (CH 3 OH, C 2 H 6 ) 
molecules through H-atom addition (i.e., reduction) reactions on surfaces of 
interstellar grains. Vertical (horizontal) dotted lines indicate “normal” abundances 
(median value for CO) as discussed in the text. For each point, ±1 a confidence limits 
are shown, and upward arrows indicate that only upper limits (3c) were obtained 
for C 2 H 2 (the arrow to the left on point 10 reflects the upper limit for CH 3 OH 
obtained for Comet S4 LINEAR). Of these comets, C/2009 Garradd (point 16) shows 
the highest CO/C 2 H 2 obtained to date, yet both its C 2 H 6 and CH 3 OH are consistent 
with their normal values. 

5-week period encompassing those observations and the ones on 
October 13 reported here. 

5. Taxonomy of primary volatile abundances: placing C/2009 PI 
in context 



I4 ‘ ‘ ‘ ‘ ‘ ‘ 0 * * * 2 * * * 4 



As a working hypothesis, comets can be classified “organics 
normal,” "organics depleted,” or “organics enriched," based on 
their abundances of the four organic molecules C 2 H 6 , CH 3 OH, 
HCN, and C 2 H 2 relative to H 2 0, with "normal” values being approx- 
imately 0.6%, 2%, 0.2%, and 0.2%, respectively (e.g., see Mumma and 
Charnley, 2011, and references therein). In C/2009 PI, the abun- 
dance ratios (relative to H 2 0) we measured for C 2 H 6 (0.64%), 
CH 3 OH (2.1%), and HCN (0.24%) were consistent with “organics- 
normal,” however C 2 H 2 (0.06%) was well below its normal value. 

The abundance of NH 3 (0.5%) fell within the range measured in 
comets (Magee-Sauer et al., 2008a), and H 2 CO (0.11%) was similar 
to that found in several comets (0.1 -0.3%); e.g., C/2004 Q2 (Mach- 
holz; Bonev et al., 2009), C/2001 WM 3 (LINEAR; Radeva et al., 
2010), and C/2007 N1 (Lulin; Gibb et al., 2012). It was however 
well below that found in C/2002 T7 (LINEAR; DiSanti et al., 2006) 
or C/2002 Cl (Ikeya-Zhang; DiSanti et al., 2002). Whether a 
“dichotomy” in (native) H 2 CO/H 2 0 (i.e., >0.5%, versus significantly 
below 0.5%) extends to the overall comet population is an open 
question that can only be answered by measuring its abundance 
in the nuclei of a significantly larger number of comets. 

Our abundance ratio CH 4 /H 2 0 in C/2009 PI (0.8-0.9%) was close 
to its median value found among comets (0.9-1%). Our measured 
abundance ratio for CO (9.1%) was within the range found in 
comets, which span <0.5% (Weaver et al., 2011) to ~30% (e.g., see 
Bockelee-Morvan et al„ 2005), but significantly above its median 
value among comets (~4%). Only C/1996 B2 (Hyakutake; DiSanti 
et al., 2003), C/1995 01 (Hale-Bopp; Bockelee-Morvan et al., 
2009; DiSanti et al„ 2001), and C/1999 T1 (McNaught-Hartley; 
Mumma et al., 2001c) had CO/H 2 0 exceeding 10%. The abundance 
of CO in C/2009 PI was also well above that found for the majority 
of comets observed within R h = 2.5 AU in the Akari survey (Ootsubo 
et al., 2012). This classifies C/2009 PI as a CO-rich comet, as noted 
from other studies (Paganini et al., 2012; Villanueva et al„ 2012c; 
Biver et al., 2012; McKay et al., 2012, Feaga et al., 2012). 



Projected Distance from the Nucleus (1 0 3 km) 


Fig. Al. Graphical illustration showing the methodology used to estimate the 
fraction of H 2 0 released in the coma of C/2009 PI (Garradd) on UT 201 1 October 13. 
We take the C 2 H 6 + HCN profile as our representative (purely) native profile (see 
Section A.1). 


6. Redox chemistry of grain surface reactions 

6.1. Molecular abundances in comets as a possible test for redox 
chemistry 

The combined reactions of condensed interstellar molecules 
(e.g., CO and C 2 H 2 ) on grain surfaces with hydrogen (reduction 
reactions) and oxygen (or OH; oxidation reactions) are known to 
produce more complex molecular species that can then be stored 
as ices in cometary nuclei. In the infrared, comparing measured 
abundances (among comets) of CO and C 2 H 2 to species produced 
through their reduction on grain surfaces (CH 3 OH and C 2 H 6 , 
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Table A1 

Estimates of water fraction released in the coma on UT 2011 October 13. 


Excess/total 

0 (excess/total) 

Mean KL1, KL2 

crMean KL1, KL2 

( l ) Summing over residual and total profiles* 


KL1 0.248 

0.0600 

0.296 

0.0296 

KL2 0.311 

0.0340 



KL12 0.289 

0.0404 



(2) Based on Q-curve analysis'* 



KL1 0.263 

0.0528 

0.284 

0.0368 

KL2 0.304 

0.0514 




a Fraction of total H 2 0 attributable to the excess emission at projected distance 
1300-1500 km sunward of the nucleus, along with its ler error, based on summing 
over its central 3000 km, corresponding approximately to the range in projected 
distance -100 to +2900 km in Fig. Al. Weighted mean from KL1 and KL2 is listed to 
the right. The entries for KL12 correspond to summing over the residual excess 
profile in Fig. Ale, representing the weighted mean of those from KL1 (Fig. Ala) and 
KL2 (Fig. Alb). The corresponding weighted mean fraction obtained by summing 
over the central 4000 km is 0.268 ± 0.049, consistent with the “Mean KL1 , KL2” 
value (0.289) listed in the table. 

b Fraction of total H 2 0 from symmetrized Q-curves. Localized isotropic release is 
assumed, as would be expected from large and therefore slow-moving particles, 
although this could be realized even for a distribution of particle sizes (see 
Appendix A). 

respectively) provides a potential means of testing the relative 
influence of oxidation. 

However, compositional measurements provide only the net 
processing history of ices as opposed to their initial abundances 
on grain surfaces. Their interpretation is therefore subject to cave- 
ats, for example: (1) Oxidation of CO leads to C0 2 , which has re- 
cently been measured with surprisingly high abundance in a 
significant fraction of comets, frequently higher than the CO abun- 
dance (Ootsubo et al., 2012). Directly measuring C0 2 requires 
observations from space and, including C/2009 PI (Feaga et al., 
2012), C0 2 has been measured in only three comets for which 
the suite of primary volatile abundances have also been character- 
ized using IR spectroscopy, the others being 9P/Tempel 1 (Mumma 
et al., 2005; Feaga et al., 2007) and 103P/Hartley 2 (Dello Russo 
et al„ 2011; Mumma et al., 2011; Weaver et al., 2011, A’Hearn, 
2011). (2) H-atom addition to CO produces H 2 CO in addition to 
CH 3 OH. Native, monomeric formaldehyde (H 2 CO) has been mea- 
sured in a number of comets, 5 however the amount contained in 
polymers is typically unknown, and where measured - in 1 P/Halley 
(Eberhardt, 1999, and references therein) and C/1995 01 Hale-Bopp 
(Wink et al., 1999) - its abundance dominates that of monomeric 
H 2 CO in the nucleus. 

A complete assessment must investigate all viable formation 
pathways for oxidized carbon compounds, including gas-phase 
production of species like C0 2 and H 2 CO (Walsh et al., 2012) and 
also addition reactions on the surfaces of (interstellar) grains 
(Oba et al., 2010; Watanabe et al., 2004). Here we discuss the po- 
tential usefulness and limitations of measured abundances in com- 
ets in the context of grain-surface reactions. 

6.2. Reduction on grains 

Hydrogen atom addition to C 2 H 2 on the surfaces of interstellar 
grains as a means of producing C 2 H 4 and then C 2 H 6 was suggested 
by Tielens (1992), and was later borne out by the discovery of 
abundant ethane in comet C/1996 B2 (Hyakutake) (Mumma 
et al., 1996). Measurement of C 2 H 6 in nearly all subsequent comets 
established this molecule as an important constituent in cometary 


5 The combination of small (sub-arc-second) beam sizes of modern IR spectrom- 
eters such as N1RSPEC, together with the relatively small distance between A- and El- 
beam positions (Section 2) favors measurement of species released directly from the 
nucleus (i.e., of “native” abundances as defined previously). 


ices (e.g., see Mumma and Charnley, 2011; DiSanti and Mumma, 
2008; Bockelee-Morvan et al., 2005, and references therein), and 
H-atom addition as both common and significant in the evolution 
of pre-cometary ices. 

Similarly, H-atom addition to CO on grains produces the highly- 
reactive formyl radical (HCO), which readily converts to H 2 CO and 
then to CH3OH. This was addressed theoretically (Charnley and 
Rodgers, 2009), and also was demonstrated observationally as a 
far more efficient means of accounting for methanol abundances 
in interstellar sources compared with gas-phase ion-neutral chem- 
istry (Wirstrom et al., 2011; Bottinelli et al„ 201 0). Laboratory irra- 
diation experiments (Watanabe et al., 2004; Hudson and Moore, 
1999) showed such surface reactions to be viable, providing signif- 
icant product yields at very low temperatures (10-20 K), but with 
virtually no production at (or above) ~25 K (Watanabe et al., 
2004). Irradiation of H 2 0 and CO ice mixtures at 10 K using ener- 
getic electrons (Bennett et al., 2011) and at 16 K with 0.8-MeV pro- 
tons (Hudson and Moore, 1999) also produced formic acid 
(HCOOH). The relative abundances of CO, H 2 CO and CH 3 OH in com- 
ets can test their inter-relationship, and the laboratory yields pro- 
vide a comparison in assessing natal conditions. 

Among the sample of comets measured at IR wavelengths, C/ 
2009 PI exhibited the most pronounced difference between abun- 
dances of CO (enriched) and C 2 H 2 (strongly depleted). Given the 
preceding discussion, this appears at odds with our results indicat- 
ing approximately “normal” abundances of C 2 H 6 and CH 3 OH. One 
possible explanation is more efficient conversion of initially less 
abundant C 2 H 2 relative to CO on its pre-cometary grains (see 
Fig. 4 of Hiraoka et al., 2000 for temperature-dependence). 

6.3. Oxidation on grains 

Perhaps a more likely explanation includes oxidation reactions 
on icy grain mantles. Surface reactions of CO with OH are efficient 
at producing C0 2 (Oba et al., 2010; Garrod and Pauly, 201 1 ). This is 
supported by the preponderance of C0 2 ice toward interstellar 
sources observed with ISO (Ehrenfreund et al., 1997; Gibb et al., 
2004), and the detection of C0 2 in 17 of 18 comets surveyed by 
AKARI (Ootsubo et al., 2012) and in three comets measured by Deep 
Impact/DIXI (9P/Tempel 1; Feaga et al., 2007, 103P/Hartley 2; 
A’Hearn, 2011, and C/2009 PI; Feaga et al., 2012). It was recently 
postulated that a significant fraction of comets from both KB and 
OC reservoirs formed in overlapping regions between CO and C0 2 
snow lines (A’Hearn et al., 2012). 

6.4. Combined reduction and oxidation ("redox") on grains 

Hydrogen-atom addition to C 2 H 2 produces the vinyl radical 
(CH 2 CH). Subjecting this to subsequent oxidation/reduction leads 
to three isomers of C 2 H 4 0: ethylene oxide (c-C 2 H 4 0), vinyl alcohol 
(CH 2 CHOH), and acetaldehyde (CH 3 CHO) (Ward and Price, 2011). 
Hydrogenation of the latter two isomers provides a viable means 
of producing ethanol (CH 3 CH 2 OH) (see Charnley, 2004). These 
studies together with our findings for C/2009 PI and perhaps other 
comets with high CO/C 2 H 2 (e.g., C/1999 T1 and C/1996 B2; see Ta- 
ble 4, Fig. 6, and discussion below) suggest that in general both oxi- 
dation and reduction should be considered. 

Our results for C/2009 PI may represent a case in point. If its 
pre-cometary ices were subjected to processing that included sig- 
nificant oxidation reactions, this could help explain its high CO/ 
C 2 H 2 abundance ratio together with its "normal” C 2 H 6 and CH 3 OH. 
The depleted acetylene in C/2009 PI could result from its conver- 
sion not only to ethane but also to the above C 2 H 4 0 isomers and 
perhaps to CH 3 CH 2 OH. Furthermore, its high C0 2 /H 2 0 abundance 
ratio (Feaga et al., 2012) suggests its “initial” CO endowment 
(i.e., the amount of CO condensed onto its pre-cometary grains) 
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may have been considerably larger than the value we observe 
(9.1%). 

6.5. Comparison among comets 

Fig. 6 compares abundance ratios CO/C 2 H 2 versus CH 3 OH/C 2 H 6 
in 16 comets identified by number in Table 4. Points well above 
the horizontal dotted line (median C0/H 2 0, “normal" C 2 H 2 /H 2 0) 
and near or to the left of the vertical dotted line (“normal” 
CH 3 0H/H 2 0 and C 2 H 6 /H 2 0) could indicate a relatively more impor- 
tant role for oxidation reactions, and correspondingly that the ices 
incorporated into such comets were subjected to processing in 
relatively oxygen-rich environments. Considering the CH 3 CHO 
observed in star forming regions (Charnley, 2004), measuring or 
stringently constraining its abundance in comets 1 having high 
CO/C 2 H 2 yet modest CH 3 OH/C 2 H 6 can test this hypothesis. However, 
it is important to bear in mind that, considering the aforementioned 
caveats (Section 6.1), comets with lower CO/C 2 H 2 must also be con- 
sidered, for example if a significant fraction of their “initial” (i.e., 
condensed interstellar) CO is incorporated into unaccounted-for 
polymers (e.g., formaldehyde polymers) and/or C0 2 . 

7. Summary 

On UT 2011 October 13 ( R h = 1.836 AU pre-perihelion) we mea- 
sured nine primary volatiles (H 2 0, CO, CH 4 , C 2 H 2 , C 2 H 6 , H 2 CO, 
CH 3 OH, HCN, and NH 3 ) in Comet C/2009 PI (Garradd). On 2012 
January 8 {R h = 1.566 AU post-perihelion) we measured three of 
these (H 2 0, HCN, and CH 4 ) and obtained meaningful upper limits 
for three others (H 2 CO, NH 3 , and C 2 H 2 ). Our measured abundance 
ratios relative to H 2 0 did not change pre- and post-perihelion, 
suggesting compositional homogeneity among species measured 
on both dates, to the depths among active regions sampled. 
However, despite its smaller heliocentric distance, Q(H 2 0) was lower 
by approximately 23% on January 8 compared with October 13, 
qualitatively similar to pre-/post-perihelion differences in water 
production revealed by SOHO/SWAN observations of C/2009 PI. 

Comparison of spatial profiles revealed differences in outgas- 
sing among volatiles in C/2009 PI. Pre-perihelion profiles for 
CH 4 , CO, C 2 H 6 and HCN showed fairly similar spatial distributions, 
suggesting these four species originated from one or more active 
regions composed predominantly of a-polar ice. Most pronounced 
was the sunward-hemisphere-facing excess in H 2 0 emission on 
October 13 having peak intensity at projected distance 1300- 
1500 km from the nucleus. This was not seen in other profiles, 
consistent with release from relatively pure, polar (water-rich) 
ice particles liberated from the nucleus, perhaps in a directed flow. 
We estimate that 25-30% of the total H 2 0 encompassed in our slit 
was released in the coma. Along with C/2007 W1 (Boattini) and 
103P/Hartley 2, C/2009 PI is the third recent comet to show evi- 
dence for (at least) two distinct phases of ice within the nucleus. 

The lack of excess H 2 0 emission together with its lower water 
production on January 8 suggest loss of the source of icy grains. 
This plus similar decreases in production of HCN and CH 4 may 
result from depletion of ice in, or covering of, one or more active 
regions on the nucleus, or perhaps more likely from pre-/post- 
perihelion differences in orientation relative to the incident solar 
radiation. 

Our results for comet C/2009 PI add a unique signature to our 
growing taxonomy of primary volatile abundances in comets. Its 
enriched CO and severely depleted C 2 H 2 , together with “normal” 
C 2 H 6 and CH 3 OH, could indicate that its ices were processed in a 


6 An abundance ratio CH 3 CH0/H 2 0 = 0.02% was reported for C/1995 01 (Hale- 
Bopp) (Crovisier et at, 2004). 


relatively oxygen-rich environment prior to their incorporation 
into the nucleus. 
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Appendix A. Distinguishing pre-perihelion sources of H 2 0 in 
C/2009 PI 


A.l. Quantitative estimate of coma source contribution 

Fig. A1 illustrates the approach we used for quantitatively esti- 
mating the fraction of water released in the coma of C/2009 PI . We 
express the H 2 0 profiles in terms of column density N be am (H 2 0 
molecules nr 2 ) averaged within a beam centered approximately 
350 km sunward of the nucleus - corresponding to the approxi- 
mate peak position of Gaussian fits to the water profiles - and hav- 
ing area (A beam ) corresponding to 2400 x 580 km at the comet. This 
is related to the “nucleus-centered" production rate (Qn C ) as (see 
DiSanti et al., 2009 for formalism): 

Nbeam =^W(^ (Al) 

^beam 

where R h is expressed in AU. 

We isolated the fraction of water released in the coma as fol- 
lows. We combined the profiles for HCN (Fig. 4b) and C 2 H 6 
(Fig. 4f), after scaling each to a common level when averaged over 
±1200 km from the nucleus. Both C 2 H 6 and HCN profiles had very 
similar growth factors (~1.5) and showed relatively large degrees 
of sunward/anti-sunward symmetry. We take this combined pro- 
file to closely approximate purely native release of primary vola- 
tiles, and refer to it as our “representative native profile.” We 
separately scale this to the total H 2 0 intensity over the range 0 
to -4000 km in the KL1 and KL2 settings; this range is indicated 
by the horizontal bar in each panel of Fig. Al. 

The residual profile obtained by subtracting this scaled repre- 
sentative native profile provides a measure of the excess H 2 0 re- 
leased in the coma and contained within our slit. We took two 
independent approaches to quantifying this fraction (excess 
water/total water): (1) We summed the residual H 2 0 profile inten- 
sities in KL1 and KL2 from approximately -100 to +2900 km, and 
compared these to the corresponding sums for the total H 2 0 pro- 
file. (2) We applied our standard “Q-curve” analysis to the residual 
H 2 0 profiles (for details of this formalism, see discussions in Dello 
Russo et al., 1998; DiSanti et al., 2001; Bonev et al., 2005), and 
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compared the terminal Q to that obtained from the total H 2 0 pro- 
file (the latter corresponding to Qt 0 t for H 2 0 in the KL1 and KL2 set- 
tings on October 13; Table 2). The resulting fractions are 
summarized in Table Al. These indicate agreement between ap- 
proaches (1) and (2), with both contributing 25-30% of the total 
water production we observe in C/2009 PI on October 13 being re- 
leased in the coma and peaking in the sunward-facing hemisphere 
at projected distance 1300-1500 km from the nucleus. Interest- 
ingly, a study of HDO in C/2009 PI with H1FI on the Herschel space 
craft provided the most consistent match to observed line profiles 
if a comparable fraction of water (30%) was adopted for release 
from icy grains in the coma (Bockelee-Morvan et al., 2012). How- 
ever, given the much larger Herschel aperture, comparisons with 
our results may not be meaningful. 

Because we must allow for possible additional release of H 2 0 in 
the coma but not encompassed in the slit, this fraction of “excess” 
H 2 0 to total H 2 0 could represent a lower bound. If larger, such 
additional release may explain the difference in water production 
reported here versus those from optical observations (Section 4.2) 
and pre-perihelion SOHO H-Lya measurements (Section 4.3). 


A.2. Possible interpretation of excess H 2 0 release 


The most straightforward explanation for the observed excess 
water in C/2009 PI is release from H 2 0-rich icy particles. These 
could arise from one or more active regions on the nucleus that 
contain an over-abundance of polar (specifically, water-rich) ice 
compared with the source(s) of activity giving rise to the ambient 
coma gas. This could be part of a directed outflow that releases pre- 
dominantly the water-rich ice component giving rise to the ob- 
served excess H 2 0. 

The excess H 2 0 peaks at a relatively localized projected distance 
from the nucleus (~1300-1500 km). If polar-ice-dominated parti- 
cles are the explanation, the distance at which they release water 
is determined by the product of their outflow speeds and their sub- 
limation lifetimes. We consider a range of particle sizes, from sub- 
micron to the largest particles that can be lifted off the nucleus by 
the sublimation process. This maximum particle size can be esti- 
mated by balancing gravitational and gas drag forces at the surface 
of the nucleus (see Delsemme and Miller, 1971): 


7mdx 


9 mv gas Z 

1671 G RcPcPg 


(A2) 


In Eq. (A2), m is the mean molecular mass of the coma gas (taken to 
be ~20 amu), Z is the gas sublimation rate (molecules s -1 cm -2 ), R c 
is the radius of the nucleus (km), and p c and p g are densities 
(gm cm -3 ) of the bulk nucleus and the particle, respectively. 
Assuming reasonable values ( R c = 2 km, a native gas production rate 
of 8 x 10 28 molecules s _1 , and p c = p g = 1 ) leads to a max ~ 5 cm. Par- 
ticles this large would remain very close to the surface of the nu- 
cleus and move very slowly. Indeed, the clumps of water-ice 
observed during the EPOXI flyby of 103P/Hartley 2 had extremely 
low velocities (<2 m s -1 ), with only one of 50 tracked in space craft 
images found beyond 15 km from the nucleus (A’Hearn, 2011). 

The smallest particles will maintain an outflow speed closer to 
(but still below) v gas , assumed to be 0.59 km s -1 at R h = 1.83 AU 
(Section 3.2), thereby requiring a lifetime of approximately 
2200-2600 s for release in the region of peak excess water emis- 
sion. (We note the relatively small phase angle of 31.6° on UT 
2011 October 13 means that if release occurs along the actual co- 
met-Sun direction, this distance and therefore the required particle 
lifetimes are approximately doubled.) 

This lifetime is consistent with grains dominated by water ice, 
as predicted by modeled sublimation of grains composed of admix- 
tures of water ice and rocky material (i.e., “dirty” grains; Beer et al., 


2006). Using a parameter X ^ (=mass wate r-ice/mass grain ), Beer et al. 
predict the lifetime of pure water ice grains (Xi = 1 ) to be about 
10 5 s at 1.8 AU, even for small grains of size a = 0.1 pm. However, 
the predicted lifetime drops precipitously when even a small frac- 
tion of refractory material is present (e.g., to 10 3 s for X , =0.9). 
Large grains (but still much smaller than a m axl e.g., a = 1 mm) hav- 
ing X! = 0.9 are predicted to survive for ~10 5 s, and these could re- 
lease H 2 0 in the coma on the observed scale provided their outflow 
speed is, as one would expect, much lower (20-30 m s -1 ). If asso- 
ciated with a nucleus rotation period of 10.4 h (see Section 3.3), the 
“distance doubling” seen in the KL2 water profile (Fig. 4a) implies a 
similar mean (projected) particle speed of ~40 m s -1 . In this man- 
ner, our observations of “excess” H 2 0 in C/2009 PI can be ex- 
plained by relatively localized release from relatively pure-water 
icy grains (0.9 <Xj < 1.0) spanning a wide range of sizes. 
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